Quantized Transport in Graphene p-n Junctions in Magnetic Field 
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Recent experimental work on locally gated graphene layers resulting in p-n junctions have revealed 
quantum Hall effect in their transport behavior. We explain the observed conductance quantization 
which is fractional in the bipolar regime and integer in the unipolar regime in terms of quantum 
Hall edge modes propagating along and across the p-n interface. In the bipolar regime the electron 
and hole modes can mix at the p-n boundary, leading to current partition and quantized shot noise 
plateaus similar to those of conductance, while in the unipolar regime transport is noiseless. These 
quantum Hall phenomena reflect the massless Dirac character of charge carriers in graphene, with 
particle-hole interplay manifest in mode mixing and noise in the bipolar regime. 



The transport properties of graphene, 2-dimensional 
sheets of graphite [l| , in particular the high carrier mo- 
bility and tunability of transport characteristics, make 
this material attractive for applications in nanoelectron- 
ics |2j, l3( • Various methods have been developed for pat- 
terning graphene sheets into prototype devices such as 
quantum dot transistors [l| and nanoribbons 0, d, fol- 
lowed by demonstration of local control of carrier density 
in a graphene sheet [6j . Besides possible device applica- 
tions, graphene junctions are predicted to host new and 
exciting phenomena reflecting massless Dirac character of 
carriers in this material, such as Klein tunneling Q, par- 
ticle collimation 0] , quasibound states 0] , and Veselago 
lensing[10]. In addition, interesting phenomena are ex- 
pected in gated graphene bilayers, where the field-effect 
transport can be induced by tuning the gap at the Dirac 
point . These applications make gating of graphene a 
topic of great interest. 

Recently, a graphene p-n junction with individual con- 
trol of carrier density in two adjacent regions with a pair 
of gates above and below it was reported [12| . The den- 
sity in each region could be varied across the neutrality 
point, allowing pn, pp and nn junctions to be formed 
at the interface. The interface width was quite small 
owing to 30 nm distance to the top gate and its sharp 
edge. Transport measurements, carried out in the quan- 
tized Hall effect (QHE) regime at fields 3 T < B < 8 T, 
revealed ohmic two-terminal conductance taking quan- 
tized values g — 6,2,3/2,1 in the units of e 2 /h. Along 
with the QHE plateaus with g = 2 and 6 in the unipolar 
regime, quantized plateaus 5=1 and 3/2 of similar qual- 
ity were observed in the bipolar regime. While conduc- 
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tance of 64- and 2^- is a hallmark of the integer QHE in 
graphene 3, 14], quantized conductance values 3/2 and 
1 are unusual and call for explanation. 

We interpret these observations by linking them to the 
properties of the Dirac-like carriers, which give rise to 
bipolar, electron and hole, QHE edge modes at the pn 
interface (FigfT]). The behavior at the interface is ex- 
plained by emplo ying ideas from the theory of quantum- 
chaotic transport [Hi EH [U El (3 0] • Although in our 
case the edge modes carry charge along the pn interface 
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FIG. 1: Schematic of QHE edge states for (A) bipolar regime, 
and (B) unipolar regime of graphene junction. In case (A) the 
edge states counter-circulate in the n and p regions, bringing 
to the pn interface electrons and holes from different reser- 
voirs. Mode mixing at the interface leads to the two-terminal 
conductance fl]). In case (B), since the edge states circulate 
in the same direction without backscattering or mixing, con- 
ductance is determined by the modes permeating the whole 
system, g — min(|i/i|, \v2\). 



all in the same direction, in a chiral rather than chaotic 
fashion, we argue that inter-mode scattering within the 
pn interface region gives rise to dynamics with features 
analogous to those known for quantum-chaotic systems. 

In this analogy the QHE states at the sample boundary 
play the role of perfect lead channels of chaotic quantum 
dots , OH , bringing charge to the pn interface and car- 
rying it away into reservoirs. However, several physical 
effects causing conductance fluctuations in chaotic dots 
are absent in our case, leading to quantization of two- 
terminal conductance not known for the dots. In par- 
ticular, the effective lead channels are quantized more 
perfectly than in the dots, owing to backscattering sup- 
pression in QHE transport. In addition, the quantum- 
mechanical interference effects which lead to sample- 
specific conductance fluctuations, can be suppressed in 
our case due to self-averaging, as well as dephasing and 
electron-electron scattering. Other effects that can affect 
the edge state transport at the pn interface are inter- 
mode relaxation and coupling to electronic states in QHE 
bulk, causing dephasing in a manner similar to the volt- 
age probe model 21 1. While these regimes yield similar 
results for conductance, they will manifest themselves 
differently in other characteristics, in particular in elec- 
tron shot noise 1 22] , which can be used for detailed char- 
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FIG. 2: Two-terminal conductance vs. gate voltage, given 
by Eq.fj] in the bipolar case v\ > 0, v-z < , and by Eq.[2] 
in the unipolar case {vi,2 of equal sign). The boundaries of 
QHE regions are specified by vip — 0, ±4, ±8..., with the gate 
voltage dependence of given by Eq.[3] Parameters used: 
distances to the top and back gates h — 30 nm, d = 300 nm, 
magnetic length Ib = 10 nm, dielectric constant k = 3. 



acterization of transport mechanisms. 

Due to particle-hole symmetry of carriers in graphene, 
the QHE in this material occurs symmetrically about the 
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neutrality point at the densities v = ±2, ±6, ±10 

In each of these quantized states there are n = \v\ 
edge modes propagating in different directions at v > 
and v < 23, 24 1. For the bipolar case, assuming QHE 
at densities v\ > and v 2 < on either side of the 
boundary, this gives \ and \v 2 \ edge modes circulating 
in opposite directions that merge to form a multi-mode 
edge state at the pn interface (FigfTJ\.). These modes sup- 
ply to the pn interface particles from both the n and p 
reservoirs. After propagating together along the inter- 
face these particles arrive at the sample boundary where 
they are ejected into the edge modes which split up and 
return to reservoirs. 

The observed conductance quantization can be readily 
explained by assuming full mixing of these modes at the 
pn interface, such that for each particle the probability 
to be ejected into any of the N — \vi \ + \v 2 | modes equals 
Pn = l/N irrespective of its origin. The two-terminal 
conductance is then obtained by multiplying pn by the 
numbers of the modes, giving 
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where v\ >2 = ±2, ±6, ±10.... This agrees with the ob- 
served quantized values (see Fig [2$. 

The character of QHE edge transport in the unipo- 
lar regime is quite different. In this case, nn or pp, the 



FIG. 3: Shot noise Fano factor, Eq.[8] plotted vs. gate volt- 
ages for the same parameter values as in Fig[2] Noise is zero 
in the unipolar regime (pp or nn) due the absence of current 
partition at the junction interface, but finite in the bipolar 
regime due to edge mode mixing at the pn interface. 



edge modes in both regions circulate in the same direc- 
tion. As a result, some modes are coupled to both reser- 
voirs, while the others are connected to only one of the 
reservoirs (FigJTp). With backscattering suppressed by 
QHE, the conductance across the boundary is solely due 
to those edge modes that permeate the entire system, 
making contact with both reservoirs. This gives the ob- 
served nonclassical conductance values 



9xm = Spp = min(|i/i|, \u 2 \) = 2,6, 10. 



(2) 



!/i : 2 = ±2, ±6, ±10..., in agreement with the known re- 
sults for quantized conductance of constrictions between 
different QHE states 25|, 2^ . The nondissipative charac- 
ter of transport in the unipolar regime, Eq.[2j resulting 
from suppressed backscattering, can be revealed by mea- 
suring noise. In the absence of current partitioning in- 
side the sample we expect only thermal Johnson-Nyquist 
noise S = 2gk^T in this regime but no shot noise contri- 
bution (see FiglHJ). 

The conductance values given by Eqns.[2] and Q] occur 
in a particular pattern [i~2| that can be described as fol- 
lows (see Figj2j. Electron density in graphene induced 
by the back gate is n\ = («;/47re)Vb/d, where d is the 
distance to the gate, Vb is voltage on it, and k is dielec- 
tric constant. Similarly, in the locally gated region we 
have n 2 = (ii/4ne)(Vb/d, + V t /h), where h and V t are 
the distance to the top gate and voltage on it. For the 
Landau level filling factors = (hc/eB)m t 2 we find 

is 1 = (£ 2 B K/2e)V h /d, v 2 = {e%K/2e){V h /d+V t /h), (3) 



with Ib = \Jhc/eB the magnetic length. The values Vb, 
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Vt corresponding to integer QHE states are inside par- 
allelograms with the boundaries approximately given by 
2/l : 2 = 0, ±4, ±8..., as appropriate for the four-fold degen- 
erate graphene Landau levels [IH, [3] ■ The resulting con- 
ductance pattern, shown in Figf2]for realistic parameter 
values, strikingly resembles the experimental results [l2j ]. 

How is the conductance in Eq.[T] affected by quantum- 
mechanical interference effects? Random matrix theory 
(RMT) of chaotic transport predicts [H, [3] ensemble- 
averaged conductance g = n\nil(ri\ + n 2 + 1 — 2//3), 
where rii.2 is the open channel number, and /3 = 1,2,4 
for the three random matrix universality classes. In our 
QHE case, with the channel numbers 711,2 — \v\,2\ and 
(3 = 2, RMT predicts g identical to Eq.[T] Similarly, semi- 
classical description of transport in chaotic cavities [13] , 
where mixing is due to the dynamics in the cavity, yields 
conductance values close to the classical result for two 
conductors connected in series. 

To clarify the origin of the mode mixing at the pn inter- 
face we studied electron density distribution for the gate 
geometry used in Ref. [T3| . Numerical solution of Laplace 
problem for the electrostatic potential in between the 
gates revealed that the pn density step is about 40 nm 
wide, a few times larger than the magnetic length at 
B = 8 T. Comparison to the known results 27] for a 
compressible region sandwiched between incompressible 
regions then suggests the presence at the pn interface of 
additional QHE modes that, in the presence of disorder, 
can facilitate inter-channel scattering and mixing. 

In the fully coherent regime conductance would exhibit 
universal fluctuations, UCF. The magnitude of UCF pre- 
dicted for chaotic transport (see Ref. [201]) in our case de- 
pends on the channel numbers as follows: 
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Applied to the observed plateaus with (^1,^2) = 
(2, -2), (2, -6), (6, -2), Eq.H indicates that these 
plateaus would not have been discernible in a system 
with fully developed UCF. We therefore conclude that 
the observed quantization of g depends on some mecha- 
nism that suppresses UCF. For example, the suppression 
could easily be understood if Thouless energy for the 
states at the pn interface was small compared to k^,T. 
The reduced UCF would then result from averaging over 
the k^T energy interval. However, the plateaus in 12] 
remain unchanged when temperature is reduced from 4K 
to 250 mK, making such a scenario unlikely. 

The UCF suppression may signal a fundamental de- 
parture of chiral QHE dynamics from that of the ear- 
lier studied systems. However, at this point we cannot 
exclude other, more mundane explanations. In partic- 
ular, time-dependent fluctuations of system parameters 
can supercede mesoscopic fluctuations, turning the ob- 
served time-averaged quantities into ensemble-averaged 
quantities. This self-averaging could arise naturally due 



to fluctuating electric field at the pn interface induced 
by voltage noise on the gates. Another, more interest- 
ing explanation could be that UCF suppression indicates 
presence of dephasing due to the coupling of the chiral 
modes to the localized states in the bulk, or some other 
intrinsic mechanism. 

Current partition due to mode mixing at the pn inter- 
face will manifest itself in the finite shot noise intensity. 
To evaluate noise, we note that mixing of the reservoir 
distributions, no matter of what origin, results in particle 
energy distribution of the form 
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which at small k^T is a double step. In anal ogy with 
diffusive systems [28| , and chaotic cavities [It], EH , this 
distribution serves as a Kogan-Shulman-like extraneous 
source of current fluctuations, 



J= I n(e)(l — n(e))de 



Wl 1 1^2 

N 2 



V f 



sd 



(6) 



We relate the noise source J to the fluctuations of the 
two-terminal current by noting that, since fluctuating 
current of intensity J is injected into each open chan- 
nel, the current fluctuations flowing into the n and p 
regions will be J± = \v\\J and J2 — \v2\J- Converting 
these fluctuations into voltage fluctuations and adding 
the contributions of the n and p regions, we find the 
voltage fluctuations induced between the reservoirs: 
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Current noise can now be obtained as S = g 2 SV 2 , where 
g is the conductance ([T]). It is convenient to character- 
ize noise by the Fano factor F = S/I, describing noise 
suppression relative to Poisson noise. We find 
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(8) 



where = 2, 6, 10.... The result © is identical in form 
to the shot noise Fano factors of chaotic cavities 17, ljj ]. 
The Fano factor values (JSJ) should be contrasted with 
F « 0.29 predicted for a p-n junction in the absence of 
magnetic field [1] . 

Another regime for noise is possible if electrons, while 
traveling along the pn interface, have enough time to 
transfer energy to each other via inelastic processes. This 
will occur if r e ; <C L/v, where r e / is the characteristic 
electron energy relaxation time, v is drift velocity and L 
is the pn interface length. (A similar regime was ana- 
lyzed for diffusive [11] and chaotic fl9j] transport.) In this 
case, the electron energy distribution is characterized by 
an effective temperature T c ff which is determined by the 
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balance of the energy supplied from reservoirs and elec- 
tron thermal energy flowing out: 



1 K \W% | 

2 H + \V2 



■K 2 d = C(2)(ki| + h|)fcl7; 2 ff , 



(9) 



C(2) = 7r 2 /6. The extraneous fluctuations, EqJSl eval- 
uated for the Fermi distribution with T = T c g, give 
J = fceTeff- Repeating the reasoning that has led to Eq.[8] 
we find the noise intensity S = gkB.T c s. This expression 
resembles the Nyquist formula, except for the factor of 
two missing because the fluctuations © occur only in 
the pn region but not in the leads. Since T e ff oc V s d, this 
noise is linear in V sc i- Similar to the T — shot noise, 
we characterize it by Fano factor F = (3F) 1 / 2 /V , with F 
given by EqJS] 

We finally note that noise can be used to test which of 
the UCF suppression mechanisms discussed above, self- 
averaging or dephasing, occur in experiment 12|. For co- 
herent transport noise exhibits mesoscopic fluctuations 
similar to UCF which can be analyzed within RMT 
framework. In the absence of time reversal symmetry, 
RMT yields ensemble-averaged Fano factor 



Wx\\V2\ 
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(see Eq. 11 in Ref. [20j). For vi, 2 = 2,4,6... this gives 
F = 4/15, 12/63, 36/143.... These values, expected when 
transport is coherent but self-averaged, are different from 
EqJ5] obtained from incoherent mixing model. 

The quantized transport observed in graphene pn junc- 
tions [l2| is of different character in the unipolar and 
bipolar regimes. In the first case, transport is dissipa- 
tionless with conductance quantized to an integer. In 
the second case, mode mixing at the pn interface cre- 
ates a situation similar to that studied in the quantum- 
chaotic transport. Conductance quantized to fractional 
values observed in Ref. 12fl then results from intrinsic or 
extrinsic suppression of UCF. These transport regimes 
can be unraveled using electron shot noise, predicted to 
be finite in the bipolar regime and zero in the unipolar 
regime, with quantized plateau structure similar to that 
of conductance. 
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